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Spherically symmetric inhomogeneous dust collapse has been studied in higher dimensional space- 
time and appearance of naked singularity has been analyzed both for non-marginal and marginally 
bound cases. It has been shown that naked singularity is possible for any arbitrary dimension in 
non-marginally bound case. For marginally bound case we have examined the radial null geodesics 
from the singularity and found that naked singularity is possible upto five dimension. 
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I. INTRODUCTION 

For the last two decades or so, gravitational collapse is an important and challenging issue in 
Einstein gravity, particularly after the formation of famous singularity theorems [1] and Cosmic 
Censorship Conjecture (CCC) [2], Also it is interesting to know the final outcome of gravitational 
collapse [3] in the background of general relativity from the perspective of black hole physics 
as well as its astrophysical implications. The singularity theorems provide us only about the 
generic property of space times in classical general relativity but these theorems can not predict 
the detailed features of the singularities such as their visibility to an external observer as well 
as their strength. The CCC, on the other hand is incomplete [4,5] as it stands because there 
is no formal proof of it in one hand and on the other there are counter examples of it. How¬ 
ever, the nature of the central shell focusing singularity depends on the choice of the initial data [6]. 

It has recently been pointed out by Joshi et al [7] that the physical feature which is responsible 
for the formation of naked singularity is nothing but the presence of shear. It is the shear 
developing in the gravitational collapse, which delays the formation of the apparent horizon so 
that the communication is possible from the very strong gravity region to observers situated outside. 

The objective of this paper is to fully investigate the situation in the background of higher 
dimensional space-time with both non-marginally and marginally bound collapse [8-10]. As it 
is only to be expected, in one way or another, these works all deal with propagation of null 
geodesics in the space-time of a collapsing dust [12]. In this context we mention that Ghosh and 
Beesham [10] have also studied dust collapse for (n+2) dimensional Tolman-Bondi space-time for 
marginally bound case (/ = 0), considering the self-similar solutions. They have concluded that 
higher dimensions are favourable for black holes rather than naked singularities. Also recently, 
Ghosh and Banerjee [11] have considered non-marginal case (/ ^ 0) for dust collapse in 5D 
Tolman-Bondi model and have shown that the degree of inhomogeneity of the collapsing matter 
is necessary to form a naked singularity. 

We show in the present paper that in non-marginally bound collapse (i.e., for f ^ 0), the naked 
singularity may appear in any dimensional space-time, but for marginally bound collapse the 
naked singularity may appear only when the space-time has dimensions upto five. 


II. NON-MARGINALLY BOUND CASE 

The metric ansatz for n dimensional space time is 
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ds 2 = - dt 2 + R dr 2 + R 2 dQl_ 2 (1) 

1 + f{r) 

with cifl 2 _ 2 = dOl + <5in 2 0i[(i0| + Sin 2 92(d9 2 +.+ 6 , m 2 0 n _3d0 2 _ 2 )] 

We consider here a dust collapse with energy-momentum tensor: 


Ty, v = p(t, r)u^u v , 


( 2 ) 


with as the n-velocity. 

Now from the Einstein’s field equations for the metric (1) with the above energy-momentum 
tensor, one can obtain (choosing 8ttG = c = 1) 

= p> 


and 


p u r) = (n-2)F'(r) 
p y > 1 2R n ~ 2 R' 


( 4 ) 


where F(r) is an arbitrary functions of r, arising from the integration with respect to the proper 
time t. 


Suppose, the collapse develops from an initial surface t = ti and the above model is characterized 
by the initial density Pi(r) = p(ti,r) and f(r), which describes the initial velocities of collapsing 
matter shells. We choose the scaling of the scale factor R such that 


R(ti,r) = r, 


( 5 ) 


so that 


Pi(r) = p(U,r ) 


n — 2 
2 


r 2 ~ n F'(r) 


Now the curve t = t s (r) defines the shell-focusing singularity and is characterized by 


( 6 ) 


R(t s (r),r)= 0 (7) 

Further, within the collapsing cloud, the trapped surfaces will be formed due to unbounded 
growth of the density and these trapped surfaces are characterized by the outgoing null geodesics. 
In fact, the apparent horizon, which is the boundary of the trapped surface has the equation 
t = t a h(f) and the scale factor at the apparent horizon satisfies 


R(tah(r),r) = [F(r)]”~* (8) 

To characterize the nature of the singularity we shall discuss the two possibilities namely, (i) 
t a h < t s (0) and (ii) t a h > t s (0). The first case may correspond to formation of black hole while 
the second one may lead to naked singularity. If the apparent horizon will form earlier than the 
instant of the formation singularity, then the event horizon can fully cover the strong gravity 
region and also the singularity. As a result, no light signal from the singularity can reach to any 
outside observer and the singularity is totally hidden within a black hole. On the other hand, in 
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the second case the trapped surfaces will form much later during the evolution of the collapse and 
it is possible to have a communication between the singularity and external observers. 

Integrating once, we have from equation (3) 


t — ti = 


(n - 1 )VF 


n _ i 2 fr n ~ 3 n -1 1 f R n ~ ^ 

r^V 2 -Fi[-, a,a + 1, — ' ] - R 1 ^ 2 F 1 [-,a,a+1,-- 


F 


(9) 


where we have used the initial condition (5) and 2 Fi is the usual hypergeometric function with 


a = | Using equations (7) and (8) separately in equation (9) we have 
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t s {r) -ti = 


(n — 1)VF 


n-1 1 fr 

r 2 2 F 1 [-,a,a+l,- — 


n —3 


( 10 ) 


and 


t a h{r ) - U 


2r " 2 1 

(n- 1)-/F 


2 Fi[-,o,a+ 


jy.n—3 

F 


2F^ 
n — 1 



a,a+ 1 ,-/] 


( 11 ) 


We note that in case of homogeneous dust, the collapse is simultaneous (Openheimer-Snyder) but 
in the present inhomogeneous model the collapse is not simultaneous (in comoving co-ordinates) 
but rather the singularity is described by a curve with starting point (to,0), which is given by 
equation (10) as 


to 


t s ( 0) — ti + 


lim 


2r~ 


(n-l)VF 


2Fi[b,a,a+ 1 , —] 


( 12 ) 


Now in order to have a finite value of the above limit, we assume F(r) and f(r) [13,14] to be in 
the following polynomial form near the central singularity (r = 0) 


F(r) = F 0 r n ~ l + F ir n + F 2 r n+l + . (13) 

and 

fir) = f 0 r 2 + fir 3 + f 2 r A + . (14) 


Then equation (5) suggests that initial density profile is also smooth at the centre and we have 


Pi{r) = po + pir + p 2 r 2 + . (15) 

with pj = (»+J-i)("-~ 2 ) j?., j = o, 1, 2,... . 

As the density gradient is negative and falls off rapidly to zero near the centre so we must have 
pi = 0 and p 2 < 0 and consequently Fi = 0 and F 2 < 0. Thus the above limit (see eq.(12)) 
simplifies to 


to — ti + -— 

(n 



2 Fi[^,a,a+l,-^-] 
z r o 


Therefore, equation (11) and (16) give (after using (13) and (14) in (11)) 


(16) 


tah to — 


2 /i 


n — 1 


( 3 » -OFF ‘ 2 


2 Fi [~, fl -f 1, a T 2, — ^-]r + 
ro 


Fo 


(n - 1 )U 0 3/2 
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2 -Fi[^, a,a+ 1,-^- 


F ° 4(n-3)(3n-7)/ o F 0 5/2 (/o + J F o ) 


(7-3n)/ 0 2 F 0 2 2J Fi[i a,a+l,-^-] 
Z 1 0 


+{-2/ 0 ^o((5 - 2n)/ 2 + 2(n - 3 )f 2 F 0 ) - 4(n - 3)f 2 F 0 (f 2 - F 2 ) + 4(n - 3)/ 0 3 F 2 


+(3n - 7)/ 2 F 2 } 2 Fi[| a + 1, a + 2, - A] j 


r 2 + 0(r 3 ) 


From the above relation we note that near r = 0, t a h > or < 0 according as 


(17) 


fi 2-Fi[^,a+ l,a + 2,-^-] < or > 0. 

Though the hypergeometric function has arguments depending on n yet 2 Fi is always positive. 
So the above restriction simply implies fi < or > 0. Now according to Joshi etal, if a comoving 
observer (at fixed r) does not encounter any trapped surfaces until the time of singularity 
formation then it is possible to visualize the singularity otherwise the singularity is covered by 
trapped surfaces, leading to a black hole. Thus /i < 0 leads to naked singularity while f\ > 0 is 
the restriction for black hole. However, the above conditions are not sufficient, particularly for 
non-smooth initial density profile. Therefore, we may have locally naked singularity or black hole 
depending on the sign of f±. 


Further, to study the effect of shear on the formation of trapped surface we first evaluate the 
shear explicitly. For n-dimensional spherically symmetric dust metric the shear scalar is estimated 

by [7] 


a = 


n — 2 ( R! 

2(n — 1) \If 



Using equation(3) we get in turn 


a = 


. - 2 [{RF' - (n - 1 )R’F} + R n ~ 3 {Rf - 2 R'f)\ 


8 (n — 1) 


(F + fR n ~ 3) 1/2 


(18) 


Since at the initial hypersurface (t = U) we have chosen R(ti, r) = r , so the initial shear tr^ is of 
the form 


— 2 \{rF' — (n — 1)F} + r" 3 (?’/' —2/)] 


8(n — 1) 


r 2 (F Jr" -3 ) 1 / 2 


Thus using the power series expansions (13) and (14) for F(r) and /(r) in the above expression 
for <x; we have 


Pi = 



2 _ fir + m (/m + F m )r m 

1 ) y/(fo + F 0 ) + fir + E™=2 m (fm + F m )r 


m 


n — 2 


8(n-l) VTo+TU L' 


fir+\2(f 2 + F 2 )~ 


f? 


2(/o + F 0 ) 


r 2 + 0(r 3 ) 


(19) 


We note that the initial shear vanishes when fi = 0, (F 2 + f 2 ) = {F 3 + f 3 ) = .=0 and hence 

even if the initial shear is zero the dust distribution may be inhomogeneous. Thus from equation 
(10) t s a function of the comoving radial co-ordinate r, so that the shell focusing singularity 
appears at different r at different instants. Also the above expression for cr, reveals that the 
existence of the naked singularity is not directly related to the non-vanishing of the shear as it 
does in the marginally bound case (see ref. [14]). 
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III. MARGINALLY BOUND CASE 

In this case (i.e., f = 0) equation (3) can be integrated out easily to give 


r " 2 1 = r V - H—1^/F(r)(t - U) (20) 

where we have used the initial condition (5). 

Suppose the radius of the spherical shell R shrinks to zero at the time t c (r) then from (20) we 
have 


t c {r) - U = 


2 r 2 


n - 1 \fW) 


( 21 ) 


Now the Kretchmann scalar 


p' 2 p p' p 2 

K = [(n - 2)(n - 3) +1] R2n _ 4Rl2 - 2(n - 2) 2 (n - 3) R2n _ 3R/2 + (n - l)(n ■- 2) 2 (n - 3)(22) 

diverges at t = f c (r) i.e., I? = 0. Thus it represents the formation of a curvature singularity at 
r. In fact the central singularity (i.e., r = 0) forms at the time 


to — ti + 


2(n — 2) 
(n - 1 )p 0 


The Kretclnnann scalar also diverges at this central singularity. 


(23) 


Now if we use the expansion (13) for F(r) in equation (21) then near r = 0, the singularity curve 
can be approximately written as (using (23)) 


t c {r) = t 0 


F n 


(n-l)E 0 3/2 


(24) 


where m > 2 and F m is the first non-vanishing term beyond F 0 . Thus t c (r) > t 0 as F m < 0 for 
any m > 2. 


To examine whether the singularity at t = to, r = 0 is naked or not, we investigate whether there 
exist one or more outgoing null geodesics which terminate in the past at the central singularity. 
In particular, we will concentrate to radial null geodesics only. 


Let us start with the assumption that it is possible to have one or more such geodesics and we 
choose the form of the geodesics (near r = 0) as 


t = t 0 + ar a , 


(25) 


to leading order in t-r plane with a > 0, a > 0. Now for t in the geodesic (25) should be less 
than t c (r) in (24) for visibility of the naked singularity so on comparison we have 


a > to and 


a < — 


Frr, 


(n — 1 )Eq /2 ' 


(26) 


Also from the metric form (1), an outgoing null geodesic must satisfy 
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dt 

dr 


= R' 


(27) 


But near r = 0, the solution (20) for R simplifies to 


R = r 




(28) 


Thus combining (25) and (28) in equation (27) we get 


aar“ -1 = 


1 - (t 0 + ar a ) - {2m +^ Fm r m (t 0 + ar a ) 


1 - 2 = 

2 


l VFo (l + ) (to + ar“) 


n —3 
n— 1 


(29) 


Now if there exists a self consistent solution of this equation then it is possible to have at 
least one outgoing radial null geodesic that had started at the singularity i.e., the singularity is 
naked. In order to simplify the above equation we shall use the restrictions in equation (26) in 
the following two ways: 


(?) a > m : 


The equation (29) becomes (in leading order) 


which implies 


aar a 1 


1 + 


2 m \ ( F m \ n ~ 1 

n-lj( 2F 0 ) 


(x — 1 T 


2 m 
n—1 


and a = 



2/(n —1) 


(30) 


(31) 


Thus for a > m we have m < -—and n > 3. 

n—o 

For n = 4, m may take values 1 and 2 and we have p\ = Fi = 0 and p 2 < 0,F 2 < 0. As a 
result, a is real and positive from equations (26) and (31). Moreover, these restrictions are already 
assumed in the power series expansion for pi(r) so that the initial density gradient is negative and 
falls off rapidly near the centre. But for n > 4, m = 1 is the only possible solution for which no 
real positive solution of ‘a’ is permissible from equations (26) and (31). Hence with the restriction 
a > m, we have a consistent solution of equation (29) only for four dimensional space-time i.e., it 
is possible to have (at least) null geodesics terminate in the past at the singularity only for four 
dimension and we can have naked singularity for n = 4. 


(ii) a = m : 

In this case equation (29) simplifies to 


mar m 


l 


F m 

'2F 0 


3 — n 



(2 to + n - 1) F m 
2 (n - 1) 



a 


2m 
r n-1 


(32) 


A comparative study of equal powers of r shows that m = ^ and a depends on F m and Fo. 

Here for n = 4, to = 3 and this situation is already discussed by Singh et al [12]. For n = 5, we 
have m = 2 and from (32) we get 
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2b 2 (Ab + £) + (26 + £) 2 = 0 (33) 

where b = -j== and f = ■&. 

V fq 

We note that for real b , we must have b < — | i.e., a <--, which is essentially the 

restriction in (26). It can be shown that the above cubic equation has at least one positive real 
root if £ < -(11 + 5\/5). Thus, if F 2 < -(11 + 5\/5)F 0 2 we have at least one real positive solution 
for a which is consistent with equation (29) (or (32)). Further for n > 5, we can not have any 
integral (positive) solution for m and hence equation (32) is not consistent for n > 5. So, it is 
possible for (at least) radial null geodesics which initiate from the singularity and reach to an 
external observer without get prevented by any trapped surface for n < 5. Therefore, naked 
singularity is possible only for four and five dimensions and for higher dimensions (n > 6) all 
singularities are covered by trapped surfaces leading to black hole. 

Further, to examine whether it is possible to have an entire family of geodesics those have started 
at the singularity, let us consider geodesics correct to one order beyond equation (25) i.e., of the 
form 


t = t 0 + ar a + dr a+0 

where as before a, d, a and /3 are positive constants. Thus equation (29) is modified to 


(34) 


aar a L +(a+f3)dr' 


a+/3— 1 _ 


1 - ^VK (t 0 + ar a + dr a+ P) - (2m +^ )f ™ r m (f 0 + ar a + dr a +P) 


i _ n—1 

1 2 


VF 0 (l + (t 0 + ar a + dr a +P) 


So for a > m (retaining terms upto second order) we have 


(35) 


aar a ~ l + (a + /3)dr“ +/J_1 = 1 + 


2 m 
n — 1 


F m 

2Fq 


2m (n —3)m 

r — + Dr a —(36) 


As before, we have the values of a and a in equation (31) and 


with 


(n — 3 )m 
n — 1 


and 


d = 


D 


( 5 —n)m 
n— 1 




3) 



— 71 + 1 


(37) 


As we have similar conclusion as before for a > m so we now consider the case a = m. But 
if we restrict ourselves to the five dimensional case then m = a = 2 and we have the same cubic 
equation (33) for b. Now /3 can be evaluated from the equation 


2 + /3 


2 

(-£-46)* 


(38) 


and we must have (3 > 0, otherwise the geodesics will not lie in the real space-time. As there 
is no restriction on d so it is totally arbitrary. This implies that there exists an entire family 
of outgoing null geodesics terminated in the past at the singularity for four and five dimensions only. 



IV. DISCUSSION AND CONCLUDING REMARKS 


In this paper, we have studied spherical dust collapse in an arbitrary n dimensional space-time. 
We have considered both non-marginal and marginally bound cases separately in sections II and 
III respectively. For non-marginal case we have seen that naked singularity may be possible for 
all dimensions (n > 4). However, to get a definite conclusion about naked singularity we should 
study the geodesic equations as it has been done in section III for marginally bound case. But 
we can not proceed further due to the presence of the complicated hypergeometric functions. 
Therefore, no definite conclusion is possible for non-marginal case. 

On the other hand, for marginally bound case we have definitely concluded that naked 
singularity is possible only for n < 5 by studying the existence of radial null geodesic through 
the singularity. This result also supports our earlier speculation (see ref. [14]). Here, we should 
mention that this result depends sensitively on the choice of the initial conditions. In fact, if we 
do not assume the initial density to have an extremum value at the centre (i.e., p\ ^ 0) the naked 
singularity will be possible in all dimensions. 

Finally, we should mention that the naked singularity described above is only a local feature, it 
is not at all a global aspect i.e., it violates the strong form of CCC. For future work it will be nice 
to consider the non-marginal case more extensively so that some definite conclusion can be drawn 
regarding the state of the singularity. Also it will be interesting to study gravitational collapse for 
perfect fluid model. 


Acknowledgement: 

The authors are thankful to the members of Relativity and Cosmology Research Centre, 
Department of Physics, Jadavpur University for helpful discussion. Part of this work has been 
done during a visit to IUCAA, Pune. One of the authors (U.D) is thankful to CSIR (Govt, of 
India) for awarding a Junior Research Fellowship. 

References: 

[1] S.W. Hawking and G.F.R. Ellis, The large scale structure of space-time (Cambridge. 
Univ. Press, Cambridge, England, 1973). 

[2] R. Penrose, Riv. Nuovo Cimento 1 252 (1969); in General Relativity, an Einstein Centenary 
Volume, edited by S.W. Hawking and W. Israel (Camb. Univ. Press, Cambridge, 1979). 

[3] For recent reviews, see, e.g. P.S. Joshi, Pramanci 55 529 (2000); C. Gundlach, Living Rev. Rel. 
2 4 (1999); A. Krolak, Prog. Theo. Phys. Suppl. 136 45 (1999); R. Penrose, in Black holes and 
relativistic stars, ed. R. M. Wald (Univ. of Chicago Press, 1998); T.P.Singh, gr-qc/ 9805066. 

[4] P.S. Joshi, Global Aspects in Gravitation and Cosmology (Oxford Univ. Press, Oxford, 1993). 

[5] C.J.S. Clarke, Class. Quantum Grav. 10 1375 (1993);T.P. Singh, J. Astrophys. Astron. 20 
221 (1999). 

[6] F.C. Mena, R. Tavakol and P.S. Joshi, Phys. Rev. D 62 044001 (2000). 

[7] P.S. Joshi, N. Dadhich and R. Maartens, Phys. Rev. D 65 101501(i?)(2002). 

[8] A. Benerjee, A. Sil and S. Chatterjee, Astrophys. J. 422 681 (1994); A. Sil and S. Chatterjee, 
Gen. Rel. Grav. 26 999 (1994); S. Chatterjee, A. Banerjee and B. Bhui, Phys. Lett. A 149 91 
(1990). 

[9] S. G. Ghosh and N. Dadhich, Phys. Rev. D 64 047501 (2001). 

[10] S. G. Ghosh and A. Beesham, Phys. Rev. D 64 124005 (2001); Class. Quantum Grav. 17 
4959 (2000). 

[11] S. G. Ghosh and A. Banerjee, Int. J. Mod. Phys. £>,(2002) (accepted), gr-qc/ 0212067 (2002). 

[12] S. Barve, T. P. Singh, C. Vaz and L. Witten, Class. Quantum Grav. 16 1727 (1999). 

[13] T. Harada, H. Iguchi and K.I. Nakao, Prog. Theor. Phys. 107 449 (2002) . 

[14] A. Banerjee, U. Debnath and S. Chakraborty, gr-qc/ 0211099 (2002) . 



